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We investigate the electroproduction of r\ mesons below a center of momentum energy of 1.6 
GeV, with particular emphasis on the roles of the iV*(1535) and iV*(1520) resonances. Using the 
effective Lagrangian approach, we show that the transverse helicity amplitude of the iV*(1535) 
can be extracted with good accuracy from the new eta electroproduction data, under reasonable 
assumptions for the strength of the longitudinal helicity amplitude. In addition, although the 
differential cross section is found to to have a small sensitivity to the JV*(1520) resonance, it is 

0^ ' shown that a recently completed double polarization experiment is very sensitive to this resonance. 
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Electroexcitation of nucleon resonances (N* states) is a clean way of studying the structures of nucleons and their 
excited states. Novel experimental facilities with polarized electron beams, available at accelerators like CEBAF at 
the Jefferson Lab (J-Lab), and corresponding developments in the polarized target technology, along with outstanding 
possibilities of the large solid-angle detectors (like the one in the Hall B of J-Lab) , make the prospects of extracting 
electromagnetic amplitudes for the N — > N* excitations as functions of four- momentum transfer, Q 2 , a realistic one. 
Indeed, the first results of such experiments are already coming on line [0J|]. 

We have explored elsewhere prospects of r\ photoproduction as a way of studying both the A r *(1535) (or Sn) 
resonance || and the difficult to access resonance 7V*(1520) B (or D13). Moving away from the real photon point, 
most of our current phenomenological knowledge about the Q 2 dependence of the Sn and D13 helicity amplitudes 
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comes from the analyses of Burkcrt || and Stoler ||. For the Sn, Q 3 A^, 2 starts to scale at about 5 GeV 2 , while for 
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0^ . the D13, the ratio 
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starts to approach unity around 3 GeV 2 , both results in agreement with the pQCD counting rules JjJ. However, the 
D13 seems to be disappearing much faster than the Sn with Q 2 . Thus, it is not clear that the D13 is, in fact, behaving 
according to pQCDl] expectations, that is, it is not certain that Q 3 ^^ is scaling at these Q 2 values. Indeed, the 
D 13 has never been studied at high Q 2 (> 3 GeV 2 ). As the Q 2 value at which the pQCD predictions become reliable 
is hotly debated in the literature ||, the exceptions to scaling should play a significant role in resolving this issue. 
As we will show, eta electroproduction can provide important constraints for the N-D13 electromagnetic transition 
amplitudes. 

Our work should also have an impact on distinguishing among the variety of baryon models that have been discussed 
in the literature, starting from the old non-relativistic model |JTo|. SU(6) analyses relativized models quark 
models involving meson and gluon exchanges between quarks ||13|| , quark models with strong emphasis on relativity 
p4|, models emphasizing large N c and so on. All these approaches are attempts to model QCD in the context 
of hadron spectroscopy. Our work here should motivate these approaches to compare with the sensitive observables 
coming from the new electromagnetic studies in the baryon sector. 

The purpose of this Letter is to investigate the role of the D13 in various eta electroproduction observables. As the 
older data jL6|, taken at Bonn, DESY, and NINA, are not of good quality, the original effective Lagrangian analysis 
jl7j did not include the D13 resonance. In addition, no D13 contribution was considered in the analysis |l| of the 
recent J-Lab data. Given that recent p/ioioproduction data |jl8f are sensitive to the D13, it is important to examine its 
role in efeciroproduction. Thus, we shall use our effective Lagrangian approach to analyze the new J-Lab data for eta 
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electroproduction at Q 2 values of 2.4 and 3.6 GeV 2 . We will determine the constraints provided by these differential 
cross section data on the D13 helicity amplitudes, Af , 2 , A^ 2 and S^ 2 , and suggest new experiments that show greater 
sensitivity to these amplitudes than the current data. In particular, based on our fits to the differential cross section, 
we predict that a recently completed double polarization experiment at J-Lab p9| is extremely sensitive to the D13 
helicity amplitudes. Finally, our work will provide a test of the model dependence in extracting the dominant helicity 
amplitude, Ai/ 2 of the Sn (A^ 2 ) from the data. 

The recent J-Lab experiment Q| measured the differential cross-section for 77 electroproduction in the W = ^/s range 
from 1490 MeV to 1590 MeV, and at Q 2 values of 2.4 and 3.6 GeV 2 . To a large extent, the data are independent of 
9 and (j>, suggesting the dominance of the Sn resonance in this reaction. However, there are minor deviations from 
angular uniformity hinting at reaction mechanisms other than the Sn contribution. Taking guidance from what has 
been learned at the real photon point [Q, we have analyzed these data using the effective Lagrangian approach, which 
consists of, in the tree approximation, the s- and u- channel nucleon, Sn and D13 exchanges, and the i-channel vector 
meson (p and to) exchanges. The nucleon, Sn and vector meson exchanges have been discussed in Ref. ||] and need 
not be discussed in detail here. The pseudoscalar 77NN coupling constant, g n NN, is taken to be 1.88, and the relevant 
combinations of vector meson couplings are taken to be 

A P < + Xu9v = 5 - 9 , 

X p g p t + \ u qX = 17.5 , (1) 

where is the V^-k coupling constant and g v and gt are the vector and tensor coupling constants to the nucleon, 
respectively. 

For the Sn contribution, we take the 77N, irN and ttttN branching ratios to be 0.5, 0.4 and 0.1, respectively, and 
the total width to be 150 MeV, which are close to the preferred values found in Ref. The PDG [^0| estimates the 
total width to be in the range from 100 to 250 MeV, and the 77N branching ratio to be between 0.35 to 0.55. In this 
work, we are primarily interested in the D13 contribution and defer discussion of the effect of these uncertainties on 
the determination of A^, 2 to a later publication. We note, however, that the quantity 

A s ,_./vZ 

(2) 
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can be determined nearly model-independently Q from the data. Thus, to a good approximation, the effect of choosing 
a different width or 77N branching ratio on A^, 2 can be determined from (g) by an appropriate scaling. Indeed, for 

the results of the three fits given in Table I of Ref. Q|, A^^tJT^ is nearly constant. 

For the D13 exchange, the strong and electromagnetic effective Lagrangians are fl7j J2l| , p2|| : 

L vNR = ^R^(Z) 75 Nd u r] + h.c. , (3) 
L 'nR = 'S-R^AYhxNF^ + h.c. , (4) 



2M 

L *nr = -^R^(X)(d x N)F» x + h.c. , (5) 

L 3 jN R = ~^R^AV)N(d x F^) + h.c. , (6) 

where M is the nucleon mass, /1 is the eta mass, and the tensor 9^(0) is defined as follows^] j25|: 

9 fl „(C)=g flu -[±(l + 2C)} 7fl7v . (7) 

The gauge couplings Gi, linearly related to the D13 helicity amplitudes, and the off-shell parameters p^| , V, X, Y, 
and Z , are a priori unknown and are determined from the fits to the data. The r/NDi^ coupling constant, gu, can be 
determined once the total width and r]N branching ratio of the D13 are specified. We take the total width to be 125 
MeV and the -qN branching ratio to be 0.1%, that is, the D13 r]N partial decay width is 0.125 MeV. According to the 



2 We take A = -1 in the spin-3/2 propagator 
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PDG |20), the D13 width is quite well determined, being in the range from 110 to 135 MeV. On the other hand, the 
77N branching ratio is quite uncertain. It is known to be small and nonzero 0, but at this point no reasonable error 
can be given to it. The 7rN and 7T7tN branching ratios are taken to be 0.6 and 0.4, respectively. 

Since the D13 is a J p =3/2~ state, the multipolarities of the jNDis transition can be electric dipole (El), mag- 
netic quadrupole (M2) and, if the photon is virtual, Coulomb dipole (CI). The relevant resonant electroproduction 
multipoles are denoted by E2-, M 2 _ and S2-, respectively, and the transverse multipoles are related to Walker's |24| 
electroproduction helicity amplitudes by, 

A 2 _ = (3M 2 _-£ 2 _)/2, 

B 2 _ = M 2 _ + E 2 - . (8) 

In terms of the electroproduction transverse helicity amplitudes, the transverse helicity amplitudes in the 'yNDia 
transition are given by, 

Af /2 = CImA 2 _ , 

^3/2 = -^CImS 2 _ , (9) 



where, 

C = 



A-KqM D Y'^ 1/2 
K C MT V 



(10) 



Here, q is the eta three-momentum in the cm frame, Mr> is the D13 mass (1520 MeV), IV the total D13 width, T v is 
the D13 — > ryN partial decay width, and K c = (W 2 — M 2 ) / (2W) is the equivalent real photon energy in the cm frame. 
All quantities are to be evaluated at W = Mo- For the Coulomb (or scalar) transition, we have 



5 ?/2 = \/^r clmS ?- > ( n ) 

again evaluated at W = Mjj. Here, k is the three-momentum of the virtual photon in the cm frame. 

The parameters of the model are gr/NN, two combinations of the vector meson couplings, the helicity amplitudes 
Ay 2 and S*, 2 , the helicity amplitudes A^, 2 , A^ 2 , and Sfj 2 , and the four off-shell parameters, for a total of 12 
parameters. The background is predominantly s- and p-wave, and controlled by g V NN, the vector meson couplings 
and the off-shell parameters. As there are tremendous correlations amongst these parameters, we hold g V NN and the 
vector meson couplings fixed at the values given above and allow the off-shell parameters to vary in the fits. As this 
is done at each Q 2 , the fitted off-shell parameters can largely compensate for any bias introduced by our choices of 
the vector meson couplings and g^NN- 

As the data have practically no angular dependence and the experiment was performed only at one e, the polarization 
of the virtual photon, no separation of A\, 2 and S[, 2 is possible. The reason for this is that in the differential cross 
section, there is no interference term between the Eq + and 5o+ multipoles^, and therefore, some arbitrary linear 
combination of these two would lead to an angular-independent differential cross section. At lower Q 2 , a transverse- 
longitudinal separation was made pq ] and it was found the the longitudinal cross section is small compared to the 
transverse one. We assume that this continues to higher Q 2 , and do fits with Sy 2 fixed either at zero or roughly^ 
10% of A s 1/2 . 

Thus, the parameters to be fitted to the data |l|] are A*, 2 , the three helicity amplitudes for the D13, and the four 

off-shell parameters, for a total of eight parameters. The results of various fits are given in Table 1 for both Q 2 = 2.4 
GeV 2 and Q 2 = 3.6 GeV 2 . We use the CERN routine MINUIT to minimize the chi-squared and the errors on the 
parameters are the so-called MINOS errors, which accounts for correlations amongst the parameters. In all cases, the 
fit to the data is excellent. In the first fit, we have fixed S^ 2 at roughly 10% of A^ 2 and allowed Ay 2 and the D13 



3 The Eo+ multipole is ~ A^ 2 i while the 5o+ multipole is ~ Sf/ 2 jlj. 

4 Specifically, since we don't expect our result for A^ 2 to be much different than found in Ref. Jl), we fix S^ 2 t° be 10% of 
Ay 2 found in that work. 
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parameters (the three helicity parameters and the off-shell parameters) to vary. The extracted A* , 2 is in excellent 

agreement with that found in At Q 2 = 2.4 GeV 2 , the extracted helicity amplitudes of the D13 all have extremely 
large errors and are consistent with zero. This is consistent with the statement made in Ref. B that the data are 
angular independent at the one-sigma level. At Q 2 = 3.6 GeV 2 , there is a slight signal for the presence of the D13. 
To examine the role of the D13 in the fit, we have turned off the D13 and refitted A^ 2 keeping S^ 2 fixed at roughly 

10 % of Ay 2 . At both Q 2 = 2.4 and 3.6 GeV 2 , the chi-squared per degree of freedom, x 2 /dgf, increases, but not 
significantly. , 2 shifts slightly upward, but within error is in agreement with that obtained from the first fit. In a 
final fit, we examined the role of Sy 2 by setting it zero and refitting Ay 2 and the D 13 parameters. The results turn 
out to be quite close to the fit with S^ 2 fixed at 10% of Ay T 

To summarize our numerical results obtained from the differential cross section, we find Ay 2 to be 50 ± 4 x 

10~ 3 GeV -1 / 2 at 2.4 GeV 2 and 35 ± 2 at 3.6 GeV 2 , in the same units. The errors here do not take into account 
uncertainties in the branching ratio, total width and mass of the Sn, which were studied in Ref. Q. Our errors are 
dominated by uncertainties in the D13 contribution, and thus should be added to the errors found in in Ref. §. Our 
results for Ay 2 are in excellent agreement with those found in Ref. At first thought, this may come as no suprise 
since we have used values for the total width and 77N branching ratio close to their preferred values. However, it should 
be emphasized that our method of analysis is quite different than the one used in [Q. In that work, the cross section 
was written as an incoherent sum of a resonance contribution and a background contribution, and it was found that 
the background contribution was less than 1% of the resonant contribution. In our work, we have a coherent sum of 
the resonance and background, and find that the background contributes at the 10% level or more. In addition, the 
role of the D13 was ignored in Ref. 

At a W of 1.54 GeV, our results are graphically^ depicted in Fig. 1 at 2.4 GeV 2 , and in Fig. 2 at 3.6 GeV 2 . The 
solid line is the result of the first fit, while the dashed-line arises when the Di 3 is turned off, everything else held 
fixed. These two lines are distinctly different and the data do seem to favor the solid line, which contains the D13. 
However, to gauge the need for the Di 3 , one should compare with the results of the second fit, that is, with the Di 3 
turned off and Ay 2 refitted to the data. This is shown by the dotted lines in Figs. 1 and 2. Looking at the graphs, 
it is difficult to tell if the data favor the solid or dotted line. At angles where most of the data exist, the dotted line 
is roughly the average of the solid line. However, at Q 2 = 2.4 GeV 2 and <f> = 90°, the solid and dotted lines are quite 
different. Therefore, data at this angle would help pin down the Di 3 parameters. 

There are two main reasons that the current data do not tightly constrain the Di 3 helicity amplitudes, in contrast 
to what happens at the real photon point pj. First, the D13 is falling faster than the Sn as a function of Q 2 , and 
thus is simply less important at these Q 2 values than at Q 2 = 0. Second, the key to pinning down the D13 helicity 
amplitudes at the real photon point is the new polarization observables in conjuction with the differential cross section 
data. At present, similar polarization observables do not exist for electroproduction. However, an experiment has 
recently been completed |0| at the J-Lab, which should have significant bearing on the extraction of the D13 helicity 
amplitudes from the eta electroproduction data. 

In the J-Lab experiment ]l9[ ], the data of which are now in the preliminary stages of analysis, both the beam and 
target were longitudinally polarized, i.e., parallel or anti-parallel to the beam direction. The polarization of the target 
was periodically flipped, resulting in the measurement of an asymmetry, which we denote as: 

a+ = l iP = l)-a(h = l,p= -1) 

et tr(h = l,p=l)+a(h = l,p=-l)' [ ' 

where h is the helicity of the incoming electron and p is the polarization of the target with z defined to be in the 
direction of the incident electron beam. Following Ref. p6[ , but correcting some mistakes originally pointed out by 
Dmitrasinovic et al. |2jj ], the differential cross section for the case of polarized beam and target can be written as 

a(6, 4>) = a + cr e + a t + a et . (13) 

The expression for the unpolarized cross section, <tq, has a standard form and will not be given here. For the other 
contributions, we work with the transverse helicity amplitudes, h±, which are related to those of Walker ]23] by 



(H 4 ±H 1 )/V2, (14) 
(#3 T H{)/y/2 . (15) 



5 The data shown in Figs. 1 and 2 represent only a small fraction of the data jfj used in our fit. 
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Walker's Hi are in turn related to the CGLN E§ F's by 



1 



Hi{6) = r sin0cos(0/2) [T 3 + T A ] 



H 2 (0) = \/2cos(0/2) 
J? 4 (0) = \/2sin(0/2) 



^2-^1 + -(1 -cos0)(^ 3 -^ 4 ) 



sin0sin(0/2) [^3-^4] 



T% + ^1 + ^(i + cos0)(.F3 + ^4) 



The longitudinal helicity amplitudes, h z Q , are related to the CGLN [E9| JPs by 



h N - - 
n — 



K 



-K 2 



(JT 7 + JT 8 ) C os(0/2) , 



-K 2 



{T 7 -T B )sin{6/2) 



Comparing with Ref. |30| , our ^1,2,3,4 are the same as their's, and our ^7,8 are related to their T^fi by, 

fco-Tr = kT§ , 
k Q T% = kF 5 , 

where ko is the energy of the virtual photon in the cm frame. 

For a polarized beam, o~ e enters the cross section, and is given by, 



a e = h-±r^2t{\ - e) sm((f>)lm{hg h* N + h F h* F ) 
K c 

= hf e 



For a polarized target, at enters and is given by 



, . P x (y/2e(l + e)em(<l>) 1 [mX 1 + e sin(20)ImX 2 

A c L 

P a (ImYi + e cos(20)Imr 2 + 2elmr 3 + y / 2e(l + e) cos(^)ImF 4 ) 



P 2 (esin(20)ImZ 2 + \/2e(l + e) sin(^)ImZi) 

Pxftx Pyfty Pzftz •> 



where 



*1 


= w 


+ Kh\ F 


x 2 


= h F h* + N 4 


- h N h* + F 
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+ h N h* F 


Y 2 


- h N h* F - 




>3 


— n, n, 




Yi 


= h$h*I - 


h F h* N 




— a ii + 


- Khf 


Z 2 


= h^hf - 


h F hf 



If both beam and target are polarized, o~ e t also enters: 



Q 

''K c 



P x W2e(l-e) cos(^)ReXi + ^1 - e 2 ReX 2 ) 



+ P y y/2e(l-e)sm((j))ReY i ) 



P z (a/1 - e 2 ReZ 2 + v / 2e(l-e) cos(<?!>)ReZi) 

h\Px fetx H - Pyfety Pzfetz) • 



(16) 
(17) 
(18) 
(19) 

(20) 
(21) 



(22) 
(23) 



(24) 



(25) 



(26) 



(27) 



The Pi in ( |25| , |27| ) are defined in a frame in which z is in the direction of the virtual photon, y is in the direction 
of k x q, and x is in the direction of y x z. Here, k is the three-momentum of the virtual photon and q is the 
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three-momentum of the eta, both in the cm frame. If the target is 100% polarized in the beam direction, then we 
find that the relevant Pi to be used in (|25| , |27| ) are 



P z = cos (3 P x = sin (3 cos (f> P y = — sin (3 sin <\> , (28) 
where (3 is the angle between the incident beam and the direction of the virtual photon; 

cos/3 =^», (29) 

with E the lab beam energy and v the energy of the virtual photon in the lab frame. For the J-Lab experiment, we 
find A^ t becomes 

Px{ftx P fetx^) ~t~ Pyifety fty) Pz{ftz "h fetz) (0(\\ 

A et — — 7 ' [ 6<J ) 

+ Je 

with the Pi given in (|2"s|). 

For the case of a pure Eq + amplitude, i.e., total Sn dominance, one obtains the simple result 



A+ = -cos(3\/l -e 2 , (31) 

which is obviously independent of 9 and 0. For the differential cross-section, it is difficult to isolate the D13 since, 
with the current statistics, a rescaling of the Eq+ can mimic (within the error bars) the effect of the D13. On the 
other hand, for A^ t a scaling of E 0+ essentially leaves this observable unchanged. Therefore, this observable should 
be able to distinguish between the fits with and without the D13. This is verified, as is shown in Fig. 3. The solid line 
is the fit with the D13, while the dashed line is this fit with the D13 turned off. The dotted line, which lies practically 
on top of the dashed line, is the best fit without the D 13 . Thus, we see that this observable is very sensitive to the 
D 13 and should provide powerful constraints on its helicity amplitudes. 

In summary, we have investigated the role of the D13 in various eta electroproduction observables, which are 
readibly measurable, for example, in Hall B at J-Lab. The signal for the D 13 contribution is very weak in the present 
differential cross section data. A full An coverage of the differential cross section should show increased sensitivity to 
the D13. Even stronger constraints on the D13 should be provided by the asymmetry At., discussed above, which has 
been recently measured at J-Lab. Regarding the Sn, we confirm the results of Ref. fjj], but the errors in that work 
should be increased slightly due to uncertainties in the D13 sector. Finally, a tranverse-longitudinal separation would 
be useful as the Q 2 dependences of the longitudinal amplitudes are also of interest as tests of QCD-inspired models, 
and of QCD itself. 

We are grateful to C. Armstrong and P. Stoler for providing us their results on eta electroproduction prior to 
publication. We also thank P. Stoler for a critical reading of this manuscript. Our research is supported by the U. S. 
Department of Energy. 
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TABLE I. The helicity amplitudes extracted from the data |l[ using our effective Lagrangain approach. At each Q 2 , the 
first row gives the results when Ay 2 and the D13 parameters are allowed to vary with S*^ fixed at roughly 10% of A^ 2 . The 
second row shows the best fit with the D13 turned off. The last row shows the best fit under the assumption S a L / 2 — 0- 
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FIG. 1. Comparison of our fits to the J-Lab data |] at W = 1.54 GeV and Q 2 = 2.4 GeV 2 . Solid line is the fit where the 
D13 parameters vary and it is assumed (S 1/ / 2 )/(A^ 2 ) ~ 10%. The dashed-line is obtained from this fit when the D13 is turned 
off, everything else held fixed. The dotted line is the best fit without the D13. 

FIG. 2. Comparison of our fits to the J-Lab data at W = 1.54 GeV and Q 2 = 3.6 GeV 2 . Curves as in Fig. 1. 

FIG. 3. Predictions for A£ t , defined in the text, at W = 1.54 GeV based on fits to the differential cross section. Curves as 
in Fig. 1. 
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